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SUMMARY 

The  problem  of  finding  the  least  favorable  configuration  for  selecting 
the  "best"  of  k populations  i.e.  the  one  with  the  largest  location  parameter 
by  use  of  six  different  two-stage  selection  procedures  is  considered.  Each 
of  the  six  procedures  consists  of  a subset  selection  (screening)  rule  at 
the  first  stage  followed  by  another  rule  ...  ed  on  (the  first  stage  and) 
additional  samples  from  the  selected  populations  to  decide  finally  which  of 
the  selected  populations  is  the  best.  In  the  indifference-zone  approach  it 
is  (or  was)  conjectured  that  the  least  favorable  parameter  configuration  is 
of  the  slippage  type.  It  is  shown  that  this  conjecture  is  true  for  four 
of  these  procedures,  for  a fifth  procedure  it  is  proved  that  at  least 
a certain  lower  bound  of  the  probability  of  a correct  selection  has  this 
property  which  is  analogous  to  the  result  of  Tamhane  and  Bechhofer  (1979) 
concerning  the  sixth  procedure. 
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1.  INTRODUCTION 

Suppose  we  are  given  k normal  populations  n(  with  different 

unknown  means  and  a common  (known  or  unknown)  variance.  If  the  experimenter's 
goal  is  to  find  that  population  having  the  largest  mean  by  using  suitably 
chosen  samples,  then  a large  variety  of  possible  sampling  plans  and  selection 
procedures  can  be  found  in  the  literature.  In  this  paper  we  are  dealing  with 
the  so-called  two-stage  procedures  of  the  following  type: 

Stage  1 : Take  k independent  samples  (X^,...,X.n)  of  size  n,  i = l,...,k, 
from  and  select  a non-empty  subset  of  these  populations 

according  to  a pre-specified  rule  S ( X ) where  X - (X^-.-.X^)  and  X.  = 
Xjj+...+X.n,  i = l,...,k.  If  the  resulting  subset  consists  of  only  one 
population,  stop  and  decide  that  this  is  the  population  with  the  largest 
mean.  Otherwise  proceed  to  Stage  2. 

Stage  2:  Take  additionally  independent  samples  of  size  m ( Y . ^ , . . . ,Y^ () 
from  those  populations  selected  in  Stage  1.  Among  the  selected 
populations  decide  finally  in  favor  of  that  population  yielding  the 
largest  Y.  (or  X.+Y^,  where  Y<  = Yii  + --+Yi,„*  1 = 1 k- 

For  convenience,  let  us  represent  the  rules  for  Stage  1 in  the  form 

S:  IR  ► ts|0  f s il,...,k)}  where  i c S means  that  the  i'th  population 

is  included  in  the  subset  of  selected  populations.  Moreover,  let  us 

represent  the  rules  for  Stage  2 in  the  form  d = lds  1 19  f s>_  { 1 , . . . ,k)l , 

2k 

where  for  every  s d^:  IR  -+  { 1 , . . . , k } and  d (j,n)  depends  only  on  those 
4^ 's  and  n^'s  with  i 6 s. 

Let  us  now  study  in  more  detail  the  four  possible  two-stage  procedures 
(S  .d,.),  a,p  * 1,2  which  we  get  after  combining  any  two  of  the  different 
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single-stage  procedures  given  below  [(S.,d^)  and  (S^,d^)  will  be  discussed 
at  the  end  of  this  section]. 

Stage  1:  For  1 t { 1 k } let 

1 c S.(X)  Iff  X.  > max  X.-c;  c > 0 fixed, 

1 ' 1 j-1 k J 

i c ( X ) iff  X.  Is  one  of  the  t largest  values  of  X1 X^; 

t t (2,. . . ,k-l  I fixed, 

[i  c Sj(X)  iff  X^  > c^;  c.t...,c^  t IR  fixed.] 

Stage  2:  For  0 f s >_  (l,...,k)  and  1 t s let 

d,  . (X,Y)  = i iff  Y . * max  Y , , 
l,S  ‘ ‘ 1 jts  J 

d,  (X ,Y ) * 1 iff  X.+Y.  - max  (X.+Y.). 

^,s  11  jts  J J 

A correct  selection  occurs  whenever  a procedure  finally  ends  up  with 
the  population  associated  with  the  largest  mean,  which  we  may  assume  to  be 
without  loss  of  generality,  since  we  obviously  are  dealing  with  permutation 
invariant  two  stage  procedures. 

To  implement  such  a procedure  one  usually  wishes  to  guarantee  that  the 
probability  of  a correct  selection  is  at  least  P*  > k"'  over  a certain  set  of 
parameter  configurations.  Now  if  the  means  say,  1 * l,...,k,  are 
restricted  to  the  condition  ^ v^-a  with  a ^ 0 fixed,  then  it 

seems  intuitively  clear  that  the  infimun  of  the  probability  of  a correct 
selection  should  occur  at  parameter  configurations  (v,\>, . . . ,v,v+a)  with 
v t IR,  which  are  called  the  least  favorable  configurations.  Since  there 
is  no  proof  for  these  conjectures  till  iow  in  the  literature  except  for  the 
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special  case  of  k = 2 populations  for  (S^.d^)  (as  we  shall  discuss  below 
more  explicitly)  we  have  tried  to  fill  this  gap  and  solve  the  problems  in 
a more  general  setup  (without  the  assumption  of  normality).  Briefly,  we 
have  been  successful  in  proving  the  conjectures  for  procedures  using 
and  but  not  for  those  using  S^. 

Discussion  of  different  two-stage  procedures: 

(Sj.dg):  (S^.d^)  has  been  studied  by  Alam  (1974)  and  Tamhane  and  Bechhofer 
(1977,  1979).  Alam  has  proved  the  conjecture  for  k = 2 and  his  subsequent 
results  are  based  on  the  assumption  that  the  conjecture  is  true  for  all  k. 
Tamhane  and  Bechhofer  (1977,  1979)  on  the  other  hand  used  lower  bounds  for 
the  probability  of  a correct  selection  which  assume  their  mfima  at  the 
desired  parameter  configurations. 

(Sj.dj)  has  not  been  studied  up  to  now.  Surprisingly,  it  turns  out 
that  it  is  even  difficult  to  prove  the  conjecture  for  this  simpler  procedure. 
Therefore,  we  propose  a lower  bound  for  the  probability  of  a correct  selection 
which  appears  to  be  quite  good  and  which  is  minimal  at  the  conjectured 
parameter  configuration. 

Remark:  Recently  Joachim  Sehr  at  the  University  of  Mainz  has  proved  that  the 
conjecture  holds  true  for  (Sj.d,,)  in  case  of  k = 3.  The  proof  is  lengthy  and 
uses  geometrical  arguments. 

(S^.d^):  (S^.d^)  and  ^.d^)  have  been  studied  by  Somerville  (1971a)  and 

(1974).  In  both  papers  he  has  claimed  that  the  corresponding  conjectures 
have  been  proved  by  Somerville  (1954),  Fairweather  (1968)  and  Somerville 
(1971b). 


Now  the  last  paper  Somerville  (1971b)  was  shown  to  be  in  error  by 
Carroll  and  Santner  (1975),  and,  in  fact,  his  method  of  proof  does  not 
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even  work  for  (S.,,dj).  Moreover,  the  "loss  function  approach"  of  Somerville 
(1954)  and  Fairweather  (1968)  Is  not  applicable  in  our  problem  since  the 
corresponding  function  W used  there  turns  out  to  be  here  an  indicator 
function  which  clearly  does  not  have  a continuous  second  derivative. 
Therefore,  the  conjectures  for  (S0,dj)  and(S^,d,)  remained  totally  unproved 
up  to  now. 

(S^,d_):  These  types  of  procedures  may  be  used  when  the  k popu  tions  are 
compared  with  a predetermined  standard  value  Vg,  say,  for  the  means.  They 
proceed  in  the  same  manner  as  the  procedures  discussed  above,  with  the  only 
difference  that  at  Staye  1 S^(X)  now  may  be  empty,  in  which  case  we  stop 
and  decide  that  no  population  is  better  than  the  standard.  The  probability 
of  this  event  is  now  desired  to  he  at  least  p*.  say,  if  vj,...,v^  - \g 
and  the  probability  of  a correct  selection  is  then  studied  over  all 

parameter  configurations  with  v^_j  * v^-.\  and  Vg  s v^.  We  will 

show  in  this  paper  that  the  infimum  occurs  at  the  point  (\>g-A Vg-.i.Vg). 

Remark:  Let  us  finally  mention  that  S^,S.,,S^  and  d are  well-established 
one-staye  multiple  decision  procedures,  studied  and  used  in  a variety  of 
papers  which  can  not  all  be  mentioned  here.  To  give  a few  references, 

Gupta  (1956,  1965)  proposed  and  studied  , Bechhpfer  (1954)  S^.  Dunnett 
(1955),  Gupta  and  Sobel  (1958)  and  Lehmann  (1961)  studied  and  Bahadur 
and  Goodman  (1952)  and  Lehmann  (1966)  investigated  d. 
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2.  GENERAL  CONSIDERAT II NS  CONCERNING  ALL  PROCIDURES 
Let  X^,  Y^,  i ■ be  indt pendent  random  variables  where  X. 

and  Y.  have  distribution  functions  t t IR  , anu  G(n-u^).  n c IR  , 

i = 1 k.  F and  G are  assumed  to  be  known  continuous  functions  and  the 

dj's  and  i^'s  represent  unknown  locition  parameter.,.  Since  we  restrict 
ourselves  to  two-stage  procedures  wl  ich  are  invariant  under  permutations 
of  the  k populations,  we  may  assume  without  loss  of  generality,  that  we 
have  6 1 1 * « « ,6|^  j ^ ^ and  1 * * * ,^|^  | "L  u k . 

Now  let  S be  any  subset  selection  rule  for  Stage  1.  Then  using  dj  or 
d,  in  Stage  2 the  probabilities  of  correct  selections  are  as  follows: 

P.($,d,)  * l P(S(X)  si / n G(n+,  , -uJdG(n)  (2.1) 

K 1 1 ,. . . ,k-l } ' IR  its  K 1 

P. (S.d  ) « V P{S(X)  • s;  X.+Y  . X.+Y  , i t si  (2.2) 

K £ *-U k-1 ) ' l i k h 

with  the  understanding  that  here  ant  in  the  sequel  s r s U { k 1 if  both  s 
and  s appear  simultaneously.  The  product  appearing  in  (2.1)  is  defined  to 
be  equal  to  one  if  s is  empty. 

In  the  sequel  let  |A|  denote  tie  size  of  any  finite  set  A.  Now  we 
state  our  main  result: 

Theorem:  For  every  A,  A ' 0,  tr  t {1,2}  and  o^,  v> q t IR  the  following  holds: 

(i)  Subject  to  e1 ok1  < ek>6  and  up...,nk_i  <.  nk-A  P^S^d,)  assumes 

its  minimal  value  at  every  parameter  configuration  (o o.i'  + a)  and 

(u, . . . ,u  ,m+a)  wl th  p , u £ IR . 

( i i ) Subject  to  the  additional  restrictions  t'k  - Sq  and  uk  > uq  P^S^.d^) 
assumes  its  minimal  value  at  every  parameter  conf i gura t i on 
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(e0-6,. . . .0^-6 ,i‘o ) and  (p , . . . ,p ,p+a)  wi th  p > m^-a. 

Proof  (first  part):  From  expression  (2.1)  it  is  clear  that  Pk ( S .d^ ) for 
every  S is  non- increasing  in  pj,...,p^  ] and  non-decreasing  in  p^.  The 

L 

same  is  seen  to  hold  true  for  P^(S,d.,)  since  for  every  (,  t IR  and 
s ( 1 , . . . ,k- 1 } PlS(X)  * s;  X..+Y.  * Xk+\’  ’ fc  *lx  = §)  also  has 
this  property. 

This  accomplishes  the  first  step  towards  a solution  of  our  problem. 

We  can  assume  from  now  on  that  pR)  = (p p,p+A)  for  some  p t IR  . 

respectively , p > p^-a  holds.  Then  (2.1)  reduces  to 


Pk(S.d.)  - l PlS(X)  * s)  / G(n+A)|SldG(„) 

s*— ( 1 k-l } ' IR 


(2.3) 


iT-1. 


= v Pik  c S(X),  | S ( X ) 1 = r)  / G(h+a)  ' 'dG(n) 
r 1 ' ‘ IR 

or,  alternatively. 


P k ( S , d ^ ) = P{k  t 


S(X)}/  ( 


G(n+A)k  'dG(n) 


(2.4) 


k-l 


+ l P(k  t S(X),  | S(X) | . r }/  G(n+A)r'1[l-G(n+A)]dG(n). 
r=  1 ’ IR 


And  (2.2)  reduces  to 


P (S,d?)  = l PtS(X)  = s;  X.+U.  < X.+A+U. , i c si  (2.5) 

K c S— 1 1 k-l)  ’ l l k k 


where  Uj are  independently  and  identically  distributed  random  variables 

with  distribution  function  G,  which  are  also  independent  of  X1>...,X(,.  (End 
of  proof's  first  part.) 

Formulae  (2.3)-(2.5)  and  the  next  lemma  will  be  used  repeatedly  in 


Sections  3 (and  4)  when  we  give  the  second  part  of  the  proof,  consisting 
of  four  versions  corresponding  to  tie  four  procedures  under  consideration. 

Lenina:  for  every  A ^ B { 1 , . . . ,k- 1 ) , r c 1 0 , 1 , . . . , | A | ) , 


dj . ,a|^|  fc  IR  and  bj 
P{  |ti  j i t A,  Xi  ' a. } j v r; 


..b|B|  t IR 

X3  ' V J ^ B) 


i s non i nc reas i ng  in  o , . 1 ..... k - 1 . 


I2.u) 


Proof:  For  r = 0 the  assertion  is  i learly  true.  For  r > 0 and  t t A, 

(2.6)  is  equal  to 

P{|{1(1  t A,  i f t,  Xi  > a^l  < r-1;  X.  b] , j fc  B,  j f i !P{X  b . I 

♦ Pi  1 1 i | i t A,i  f < , X . ^ a . } | = r ; i bj  «j  t B,  ,i  + t)P{ - min(ae  ,b{ ) I 

which  obviously  is  nonincreasing  in  d{.  Similarly  one  can  prove  the  assertion 
in  case  of  c t B \ A,  whereas  in  case  of  c it  is  trivially  true  .since  in 
that  case  (2.6)  does  not  even  depend  on  o , . 

3.  THE  SECOND  FART  OF  THE  PROOF 
3.1  Case  (S^.d^) 

For  every  fixed  t t (2,...,k-l),  we  have  |S.,(X)|  = t with  probability 
one,  and  thus  (2.3)  reduces  to 

P|c(s2.d1 ) - Plk  t S2(X»  ^6(n+A)t‘1dG(n).  (3.1) 


Moreover,  we  have 
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P{k  € S2(X)}  = P { | { i | X k < X1)|  < t-1]  (3.2) 

= J P{ | { i | X . > 4,  i f k} | < t-1 |Xk  = C>dP{X.  = 4> 

IR  1 k k 

= / Pt  | { i | X > 4+0k,  1 + k)|  < t-1)  dFU). 

IR 


Since  the  integrand  obviously  is  nondecreasing  in  ek  and  by  the 
lemma  is  nonincreasing  in  e.| , . . . ,ek_-| , the  proof  for  (S2>d-|)  is  completed. 

3.2.  Case  (S2,d2) 

Let  t £ {2,...,k-l}  be  fixed.  Then  using  the  fact  that 
Uk  -|_uk)  is  symmetrically  distributed,  from  (2.5)  we  get 

Pk(S2,d2)  = (3.3) 


l—i  1 < 1 2<  * • - < i t_ i 


P{X*  < V 

<k-l  l ’ 

- . (:S  its 


uj-uk+xi  < Xk+&,  J*1 t-1) 

J 


s {i1,....1t_1} 


/■•■/  I.  ...  I . . . I ptx,  < V *+v 

{a1  <a2< . . . <a{._1 } IR  Ui^i^. . .<it.1<.k-l  * 

s = {ir...,it_1} 

a +X.  < 4+9.+A,  j=l t-l}dF(f  )dP{U,-U.=a.,  j=l ....  ,t-l } 

^ ij  k J k j 

where  in  the  second  sum  n runs  over  all  (t-1)!  permutations  of  (1 ,2, . . . ,t-l ) . 
Since  every  probability  term  obviously  is  nondecreasing  in  ek>  it  follows 


that  P)|<(^2’c*2^  a^s0  ^as  Pr°Perty- 
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Now  we  show  that  for  every  fixtd  f,  fc  IR  and  a^a.,  >....*  a^  ^ the 
integrand  is  nonincreasing  in  • For  sake  of  simplicity  we 

prove  it  for  e^.  Let  = c+0k+A-a  , j=l,...,t-l  and  • k k . 


l<i l<i^< . 


•<it-lik_1  " 


1 P»xt  < 

4s 


i ' ; k ’ 
its 


ni Vr 


X.  < b , j = l t-lt  (3.4) 

J "j 


= r 


j^k-2  1 <_1 1 < 1 2<  * • - i t_  1— *t“  " 

s = i i 1 . . ,it_1 ) 


l PiX,  < X.,  £.; 

i ik 

e^s  its 


X1  < V*  j = 1 . ,t-l  |X2=c2 Xk-rck-l,dPlX2=t2 Xk-rf’k-l' 

J J 


Let  now  £.2,...,t.k  j t IR  be  fixed  and  assume,  without  loss  of  generality, 
that  l2  < <...<  holds.  Then  the  integrand  reduces  to 


l P(Xr  Xk_t  < xk_t+1,  ck;  xk_t+1  < b , 


Ck-t+2  - bn2 xk-l  1 b"t.JX2  = f'2*‘  ’ ’ ,Xk-l  = ' k-1 


- PlXk-t+l  < Xl’  ‘‘k’  X1  - bn* 

TT  I 


(k-t+2  - brr2 xk-l  - b"t.  = ;2 Xk- 


1 " ' k-1 


Finally  we  have  to  distinguish  between  two  cases  according  to  whether 
^k-t+1  * ^k  or  not’  *n  case  f’k-t  + l < f'k’  (3.5)  educes  to 


(3.5) 


in  whereas  the  second  actor'  does  not  depend  on  ^ 

jnd  is  nondecreasing  in  t>^.  Thus  the  proof  for'  (S^.dj)  is  completed. 

3.4.  Case  (S-^.d,,) 


For  every  fixed  ck>  froii  (?.b)  we  have 


1? 


P.  (S,,d„)  - l IMX,  > c.,  X.  v c,; 

* J C Ci  J ) l m ] ) 1 1 J J 

11  i c s j y s 


i3.d) 


Xi  * V V A + Uk» 

i c s 

which  clearly  in  nondecreasing  in  o(  since  in  every  summand  we  have 
s * s i (XI  by  our  convention.  Aga  n,  for  sake  of  simplicity,  it  will  be 
shown  (hat  i’k(S?,d,)  is  nonincreasing  in  Oj.  Now 


WV  * ^(,  ■ k , m cr  X1  ' V xj  ' Y 

^ ‘ i ^ S J ( s 

J i 1 

x<  4 Ui  ■ Xk  4 ■' 4 V 

i c s 


(3.10) 


♦ [ PiX  ' c..  X.  ' c.,  X c,; 

s>--|  2 k-1)  ’ ’ 1 1 J > J 


r c s 


X,  ♦ u,.  Xt  ♦ U,  . Xk  ♦ A ♦ Ukl 
t c s 


j i S 

J 1 1 


PU,  ' v , * X.  < it  i X,  ' i j » 
ip  i.u  1 11  1 J J 
L 9 . • . »*  I ? a * * 1 i ^ 


\ C S 


*1  4 ui  V M UK  ' C1  4 V 

i c s 


J t 1 
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+ l PlX.  > cr  X < c ; 

swi2,...,k-n  i fc  s . i s 

j t 1 


C1  ♦ Ur  X}  + Ur  X.+U.  < Xk  + A 


1 fc  s 

Clearly,  every  summand  is  nonincreasing 
(S^.d,,)  is  completed. 


+ v* 


in  9, 


and,  therefore,  the  proof  for 


3. Concluding  Remarks 

In  the  case  of  normal  populations  as  described  in  Section  1 we  have 

2 2 2 
X.  i.  N(nv.,  rut  ) and  Y.  ^ N(mv.,  mo  ),  i = l,...,k,  for  some  o >0.  Thus  for  every 

i c 1 1 ,2,31  and  p c 11,2}  the  probability  that  (S^.d,)  finally  leads  to 
a decision  in  favor  of  population  can  be  represented  by  a certain 
function  ,(v^ ,v^,. . . ,v^).  To  prove  the  conjectures  we  could,  alterna- 
tively, have  tried  to  show  that  H , is  nonincreasing  in  , and 

nondecreasing  in  v^.  But  this  turns  out  to  be  a very  difficult  and 
cumbersome  way. 

4.  A LOWtR  BOUND  FOR  (Sj.d,) 

To  prove  the  conjecture  for  (S^,d^)  in  view  of  (2.4)  it  would  suffice 
to  show  that  for  every  r c { 1 ,. . . ,k}  Pik  cS^(X),  |o^(X)|  < r ) is 
non increasing  in  j and  non-decreasing  in  t>k.  For  r 3 1 and 

i k this  probability  is  equal  to  P(Xj,...,X^  ^ < X^-cl  and  P{X^,...,X^  ^ <_ 

Xk  * cl.  respec  ti  vly,  wlu'i  o each  of  them  clearly  has  the  desired  property. 

Thus  the  conjecture  for  k - 2 is  proved. 

lor  k 2 it  turns  out  to  be  rather  difficult  to  prove  the  conjecture. 


Therefore,  we  dorh  > a lower  bound  for  the  probability  of  a correct 
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selection,  which  Assumes  its  minimal  value  at  the  desired  parameter 
configuration. 

For  k > 2 and  r t (2 k-1)  we  have 


Plk  fc  S, (X) . |S1(X) J < rl  « 


(4.1) 


= / P<k  Vxi Vt.O,  | S, (X, X(<_1  ,c ) | . r|Xk  = U 

dP{ Xk  = U 

1 | P(|{1|Xj  > C+ok-c,  i f k}|  < r-1 , 

X j • • ■ • ,Xk  l < 4 + ek+cl  dF(t). 

The  preceding  inequality  holds  also  for  r ' 1 and  follows  from  the  fact 
that  we  have 


{k  t S ( X ) , | S( X ) | v r)  « (4,2) 

* (|(i|X.  >,  Xk-c,  i f k > | < r-1,  X.  = max  X.) 

k j = l k J 

k-1 

■ U t|(i|Xi  > X -c,  i f k,U|  < r-2,  X.  > X -c,  X,  * max  X.). 
1=1  k t y j = 1 kJ 

kow  the  integrand  In  the  last  Integral  of  (4.1)  clearly  is  nondecreasing  in 
t>k  and  by  our  lemma  It  is  nonincreasing  in  Oj ,. . . ,ok_j . Thus,  using  (2.4), 
we  get 


Pk(S1.d1)  > P{Zr. 


«.  Zk+iS+c } 


^6(n+A)k  ^dG(n) 


(4.3) 


k-1 

♦ [ P( | (1 |Z,  > Z.+6-c,  1 f k>|  < r-1, 

r«i  ix 


- -----  — 


I 


I 


1 


1 


, v Z, +-6+C  }/G(n+A)r  1-G(n+A) JdG(n) 

' k k II' 

where  Zj,...,Z^  lire  independently  and  identically  distributed  random 
variables  with  distribution  function  F.  Now  the  right  hand  side  of  (4.3), 
being  in  a form  similar  to  (2.4),  can  be  brought  into  a form  similar  to 
that  ot  (.’.3).  lhen  it  is  equal  to 
k 

l PI U i | Z . > Z.+6-c,  i f k > * r-1 , (4.4) 

r*  1 

Zj,...,Z^_j  v Z^+A+c I / G(n+A)r  'dG(u) 

IK 

* l j (*"  j ' )F((.+a-c  )^  ’ ' LP(t+A+c  )-F (t+6-c)  j'dF(c )/  G(n+A ) 'dG( n)  • 
1-0  IK  IK 

Thus  we  finally  arrive  at  the  following  result: 

Corollary:  For  k > 2 

Pk(Srd])  > (4.5) 

' / / iF(f,+A-c)+LF(c+A+c)-F(f.+A-c)]G(n+A)  'dF(f.)dG(n)- 

IK  IK 

Note  that  for  k 2 this  lower  bound  for  the  probability  of  a correct 
selection  is  exact. 

This  research  was  supported  by  the  Office  of  Naval  Research  Contract 
N00014-75-C-0455  at  Purdue  University. 
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